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1.  INTRODUCTION 


The  first  of  the  DNA  series  of  IEMP  benchmark  code  calculations 
involves  the  transport  of  electrons  in  a simple  cylindrical  cavity  under  both 
vacuum  and  finite  gas  pressure  conditions.  The  purpose  of  this  problem  was 
to  check  the  various  DNA  SGEMP  codes  in  performing  calculations  when  the 
manner  in  which  they  were  performed  was  as  close  as  possible.  The  electron 
emission  data  was  specified  by  HDL.  It  was  particularly  simple  to  model 
and  was  not  based  on  experimental  data.  In  addition,  the  grid  spacings  and 
time-steps  were  specified  so  that  a close  comparison  of  the  calculations  should 
be  possible.  This  is  discussed  in  the  following  section.  The  results  of  the 
calculations  are  presented  in  Section  3 and  these  are  then  summarized. 

To  carry  out  the  calculations,  we  used  the  S:*  codes  FRED/EM  and 
GFRED.  The  first  of  these  is  a 2-D  axisymmatric  electromagnetic  particle 
code  and  is  described  elsewhere^.  It  was  used  for  the  vacuum  case  cal- 
culations. In  order  to  handle  the  relativistic  dynamics  of  the  more  energetic 

electrons  a relativistic  particle  pusher  was  included.  It  is  based  on  the 

(2) 

method  of  Boris'  . 

The  non-zero  pressure  calculations  were  performed  using  a new 
code,  GFRED.  It  is  an  adaption  of  FRED/EM,  in  which  the  secondary 
electron  dynamics  and  ionization  processes  are  included.  It  is  described  in 
the  Appendix. 
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PRECEDING  PAGE  BLANK- NOT  FILMED 


2.  PROBLEM  SPECIFICATION 


The  problem  eonfigu ration  is  shown  in  Figure  1.  It  consists  of  a 
cylindrical  tank  whose  sides  are  conducting  and  grounded  to  each  other. 
The  radius  R and  the  length  L are  both  15  cm.  The  electron  emission  is 
uniform  over  the  omitting  surface  z 0 and  contains  no  angular  spread. 
The  spectrum  of  the  electrons  is  shown  in  Table  1 and  is  uniform  in  time. 


TABLE  1 


Mean  Energy  (keV) 
15 
25 
35 
45 
GO 


Fraction  of 
Total  Charge 
0.23 
0.34 
0. 1G 
0.13 
0.14 
1.00 


The  emitted  current  pulse  has  an  amplitude 


J(t)  =J  f-  , 
° lR 


t < t 


R 


t > t 


R 


where 

J =5  A/cm" 
o 


t = G nsec  . 
R 


The  electric  and  magnetic  fields  are  calculated  at  the  positions  indicated  in 
Figure  1. 
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Figure  1.  Cavity  geometry,  showing  approximate  positions  where  the  desired 
output  quantities  are  calculated  (accurate  to  dimensions  of  a cell 
Ar  = 1 cm,  Az  = 0.  5 cm).  Exact  positions  are  indicated  on  the 
respective  output  plots. 
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Before  proceeding,  it  is  worth  considering  the  basic  phenomena  in- 
volved. In  the  zero  pressure  case  the  emitted  current  will  be  limited  by  its 
own  space  charge  il  it  exceeds  the  (non-relativistic)  Child- Langmuir  current 
limit,  namely 

„ 3/2 


J^,  2.  33  x 1 0 


A/cm’ 


where  e is  the  emitted  electron  energy  (in  eV).  This  is  only  a rough  esti- 
mate, and  ignores  electromagnetic  cavity  resonances  as  well  as  electron 
inertia  effects,  2-D  radial  spreading  and  polychromaticity  in  the  electron 

emission.  We  shall  make  one  further  assumption,  namely  that  the  limiting 

4 

current  is  determined  by  the  most  energetic  electrons.  Taking  e 0x10 
eV  and  ignoring  relativistic  corrections,  we  see  that 

J ^0.6  A/cm2  . 

This  corresponds  to  less  than  one-eighth  of  the  peak  emitted  current.  Con- 
sequently, we  should  expect  limiting  to  occur  fairly  early  in  the  pulse. 

There  is  no  very  low  energy  component  of  emitted  electrons  in  this 
problem  and,  as  we  shall  see,  limiting  does  not  occur  very  close  to  the 
emitting  surface.  The  cell  size  used  was 


Ar  1 cm 


Az  0.  5 cm 


and  this  proved  to  be  quite  adequate.  The  codes  used  were  electromagnetic 
and  the  time- step  stability  criterion  demands  that 


C^)2  > + 

cAt  2 2 

(Ar)  (Az) 


At  < 1.5x10  secs  . 


Since  some  o|  the  electrons  have  velocities  comparable  with  the  velocity  of 
light  (3  % 0.  15),  it  is  important  they  be  pushed  using  a time-step  which  is 
comparable  with  the  above.  In  FRED/EM,  the  two  time-steps  were  identical. 
In  the  original  specifications  of  this  problem,  particles  were  permitted  to  be 
pushed  with  a much  longer  time-step  ( 5x10  secs).  With  such  a time- 
step.  particles  cross  many  cells  in  a single  cycle  and  this  will  lead  to  noisy 
results. 

The  vacuum  case  is  essentially  quasi-static  since 
el  150  cm 

and  this  is  much  larger  than  the  radius  or  length  of  the  cavity.  Cavity  reson- 
ances will  be  of  small  amplitude  and  high  frequency.  However,  the  advantages 
to  be  gained  by  using  a quasi-static  code  are  not  very  great  for  this  problem 
since  there  are  many  fast  moving  electrons  and  the  time  savings  are  only 
about  /3  1 % 2.  In  any  ease,  FRED/E  was  not  used  as  the  point  of  these 
calculations  was  to  compare  c<xles  performing  under  closely  similar  running 
conditions. 

in  the  non-zero  pressures  cases  (50  and  200  mtorr),  the  production  of 
secondary  electrons  leads  to  a reduction  of  the  space-charge  potential.  Sec- 
ondaries arc  accelerated  away  from  such  regions  and  the  beam  electron 
charge  is  largely  neutralized  by'  the  remaining  ions.  The  time  taken  to 

(3) 

neutralize  the  beam  space-charge  is,  very  crudely,  given  by' 

1 

lN  n erg 
o c 

where  n^  is  the  gas  number  density,  a is  the  ionization  cross  section  and  6c 

is  the  average  electron  velocity.  At  a pressure  of  50  microns  (n  n, 

15  -3  0 

1.8  x 10  cm  ) for  40  keV  electrons 

“8  3 

<t6c  'v  2,  9 x 10  cm  ,/sec 
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so  that 


t 'v  19  nsec. 

N 

Thus,  gas  neutralization,  while  important  for  the  50  mtorr  case,  should 
dominate  the  behavior  in  the  200  mtorr  case. 

When  neutralization  occurs,  the  beam  electron  number  density  is 
similar  to  that  of  the  secondaries.  With  beam  electrons  moving  rapidly 
through  the  cold  secondary  electron  plasma,  it  is  possible  for  collective 
instabilities  to  develop. 
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3.  RESULTS 


In  this  section,  the  results  obtained  for  each  of  the  three  calculations 
are  described  in  turn. 

P ^ 0 

Figures  2 through  9 are  snapshots  at  nanosecond  intervals  of  the 
electron  dynamics.  Deceleration  of  the  beam  is  evident  at  t = 2 nsec  and 
trajectory  turnaround  is  important  by  t = 3 nsec.  Radial  spreading  of  the 
beam  is  greatest  far  from  the  axis  of  the  tank  where  the  peel-off  of  the 
lower  energy  groups  is  most  obvious.  It  is  also  clear  that  the  spatial 
zoning  is  appropriate  to  this  problem  and  that  no  limiting  occurs  within  the 
first  cells. 

Time  histories  of  a number  of  electric,  magnetic  and  current  field 
quantities  at  points  specified  in  Figure  1 are  shown  in  Figures  10  through 
20.  In  Figure  10,  the  peak  magnetic  field  at  (r  = 15  cm,  z 15  cm)  is 

B - 5 gauss 

This  corresponds  to  a total  current  I transmitted  to  the  z -15  cm  surface 
given  by 

I = 5 r B 

= 375  Amperes  . 

The  average  transmitted  current  density  on  this  surface  is  then 
j ^ 0.  5 A/cm^ 

This  result  is  close  to  the  Child- Langmuir  result  (Equation  1),  especially  in 
view  of  the  various  assumptions  made.  Probably  the  largest  error  made  in 
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Figure  15.  Radial  electric  field  at  r = 14  cm,  z = 7.25  cm. 
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Figure  18.  Axial  electric  field  at  r = 14.5  cm,  z = 14.5  cm. 
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Figure  19.  Net  electron  current  density  on  emitting  surface  at 
r = 0.5  cm. 
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this  estimate  is  the  neglect  of  the  2-D  effect  of  radial  blowup. 

Figures  19  and  20  show  the  net  axial  electron  current  on  the  emitting 
surface  z 0 at  the  positions  r = 0.5  cm,  r 14.5  cm,  respectively.  In 
the  former  case,  the  current  grows  linearly  until  current  limiting  sets  in. 

In  the  latter  case,  the  current  is  not  limited  on  account  of  2-D  effects;  namely, 
smaller  axial  fields  close  to  the  curved  edge  of  the  tank  at  r R.  In  both 
cases,  the  sharp  jumps  in  the  current  result  from  the  return  of  a discrete 
energy  group.  With  the  lack  of  angular  spread  in  the  emission,  these  groups 
are  quite  obvious.  They  would  not  show  up  in  more  realistic  emission  data. 

The  high  frequency  oscillations  in  the  current  responses  result  from  particle 
fluctuations.  Particles  are  emitted  only  every  other  cycle  and  this  gives  rise 
to  the  very  high  frequency  noise.  It  is  removed  entirely  by  emitting  every 
cycle.  In  this  calculation,  a maximum  of  6000  particles  was  permitted. 

The  initial  rise  of  the  magnetic  field  (Figures  10,  14)  is  due  to  the 
displacement  current,  not  the  electron  current.  After  limiting  has  set  in, 
the  magnetic  field  shows  some  oscillatory  behavior.  This  results  from  cavity 
mode  oscillations.  They  are  relatively  weak  because  the  dimensions  of  the 

cavity  are  small  compared  with  ctn  so  that  the  response  is  basically  quasi-static. 

K 

P = 0.  05  torr. 

A series  of  graphical  plots  which  is  similar  to  the  P = 0 case  is 
presented  here  (Figures  21-28).  As  before,  limiting  begins  at  about  t = 3 nsec. 
However,  it  is  much  less  severe  as  the  particle  plots  show.  The  space  charge 
barrier  is  partially  neutralized  and  transmission  across  the  cavity  is  more 
complete.  Time  histories  of  electric  and  magnetic  fields  are  shown  in  Figures 


29-39. 


According  to  the  estimate  given  in  Section  2,  neutralization  should  not 


be  complete  by  the  end  of  the  problem  (8  nsecs)  since  t ^19  nsecs  at  this  pressure. 
This  result  is  borne  out  by  the  results  obtained. 


31 


0 5 10  15 

Z (cm) 


Figure  21.  Particle  positions  at  t = 1 nsec,  P = 50  mtorr. 
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Figure  23.  Particle  positions  at  t = 3 nsec,  P = 50  mtorr 
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Figure  24.  Particle  positions  at  t = 4 nsec,  P = 50  mtorr 
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Figure  31.  Radial  electric  field  at  r = 14.0  cm,  z = 0.25  cm 
P = 50  mtorr. 
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Figure  34.  Radial  electric  field  at  r = 14  cm,  z = 7.25  cm, 
P = 50  mtorr. 
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Figure  39.  Axial  electron  current  density  on  emitting  surface 
at  r = i 4.5  cm,  P = 50  mtorr. 
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The  peak  transmitted  current  is  (see  Figure  29) 

I = 5 r B (gauss)  Amps 
= 900  Amps 

or  an  average  current  density  of  approximately  1.3  A/cm“.  It  is  interesting 
to  note  that  the  current  transmitted  continues  to  increase  after  the  time 
t = 6 nsec. 

P = 0. 2 torr. 

This  series  of  graphical  plots  (Figures  40-47)  are  similar  in  form 
to  those  of  the  preceding  cases.  As  before,  limiting  is  in  evidence  by  3 nsec. 
However,  it  is  even  less  severe  in  this  case.  The  position  where  turnaround 
occurs  reaches  forward  toward  the  opposite  face  as  the  gas  becomes  increas- 
ingly ionized  and  large  numbers  of  secondaries  are  driven  away  from  regions 
of  high  potential  to  reduce  the  space-charge.  Late  in  the  problem  a large 
number  of  secondaries,  which  now  dominate,  are  subject  to  drift  motion  in 
the  magnetic  field  as  quasi-neutrality  has  been  attained. 

The  total  current  transmitted  (from  Figure  48)  reaches  a large  value 
by  6 nsec 

I = 75  B 

= 2200  Amperes 

2 

or  an  average  current  density  of  over  3 A/cm  . This  amounts  to  over  half  the 
injected  current.  The  result  appears  to  be  consistent  with  those  of  the  previous 
cases. 
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Figure  41.  Particle  positions  at  t = 2 nsec,  P = 0.2  torr. 
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Figure  44.  Particle  positions  at  t = 5 nsec,  P = 0.2  torr. 
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There  are,  however,  some  features  of  the  P = 0.2  torr  which  are  quite 
different  from  the  two  other  cases.  On  account  of  the  higher  gas  density, 
quasi-neutrality  is  reached  before  the  end  of  the  calculation.  With  energetic 
beam  electrons  moving  through  a cold  plasma  of  secondaries  having  a com- 
parable number  density,  conditions  are  suitable  for  collective  beaming 
instabilities.  In  fact,  at  t = 5 nsec,  the  spatial  profiles  of  the  electric- 
fields  exhibit  large  amplitude  variations  (not  shown  here). 

Later  in  the  calculation,  after  6 nsec,  the  magnetic  field  shows  large 
fluctuations  of  high  frequency  (Figures  48  to  50)  which  continue  to  grow  in 
amplitude.  While  the  fluctuations  may  result  initially  from  real  collective 
effects,  they  are  certainly  not  correctly  treated  at  the  end  of  the  calculation. 
By  that  stage,  the  stability  criterion 


co  At  < 1/2 
P 


is  violated,  resulting  in  unstable  electrostatic  oscillations.  Here, 
electron  plasma  frequency  given  by 


co  is  the 
P 


4 1/2 

co  'v  5. 6 x 10  n 
P 


The  electron  number  density  n is  given  by 


nep  = 0. 1 j 

-10  2 

where  e = 4.8x10  esu  and  j is  in  A/cm  . In  the  vacuum  case,  we  have 
, 2 

j 'v  0. 5 A/cm 

8 —3 

and,  with  0^1/2,  we  find  that  n % 2 x 10  cm  , or 

Q -1 
O0p  * 8 x 10  sec  , 
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which  clearly  satisfies  the  stability  condition, 

_ o 

In  the  finite  gas  density  case,  however,  many  low  energy  (/j'vlO  ) 

2 

electrons  are  produced  and  the  current  flowing  is  larger,  3 \ 5 A/cm  , so  that 

,„12  -3 

n 'v  10  cm 


or 

„ ,10  -1 
u)  'v  5 x 10  sec 
P 

In  this  latter  case,  the  time-step,  at  least  late  in  the  problem,  is  required  to 
be  as  small  as,  or  even  smaller  than,  the  electromagnetic  time-step.  This 
effect  is  discussed  further  in  the  Appendix. 
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4.  SUMMARY  OF  RESULTS 


The  vacuum  case  results  were  quite  satisfactory  and  should  provide 
a useful  basis  for  comparison  with  other  code  calculations.  Even  with  the 
relativistic  particle  pusher  included  the  computation  time  of  2.5  sec/cycle 
on  the  Univac  1108  for  approximately  4000-5000  particles  is  not  unreason- 
able, namely  500  nsec  per  particle  per  cycle.  This  is  over  twice  as  long  as 
the  non  relativistic  pusher.  The  results  obtained  were  smooth  and  free  of 
any  large  statistical  fluctuations. 

The  P =0.05  torr  case  also  provided  some  interesting  and  useful 
results.  This  was  the  first  use  of  the  finite  pressure  code  GFRED.  It 
appeared  to  behave  in  a satisfactory  manner  throughout  with  computation 
times  similar  to  those  of  the  pressure  treatment  in  FRED/EM.  There  were 
no  unusual  statistical  fluctuations  and  the  results  appeared  to  be  consistent. 

The  P =0.2  torr  case  provided  a more  severe  test  of  the  GFRED 
treatment.  It  behaved  in  a satisfactory  manner  for  most  of  the  calculation. 
However,  after  rapid  multiplication  of  the  secondaries  had  begun  it  was 
evident  that  abnormal  fluctuations  were  present.  This  problem  could  have 
been  solved  by  reducing  the  time-step,  but  for  the  present  work,  this  was  not 
done.  It  appears  that  an  extended  version  of  GFRED  is  required  if  space- 
charge  neutral  cases  are  to  be  examined.  The  existing  GFRED  treatment 
works  well  when  the  electron  density  does  not  impose  severe  time  constraints 
but  for  high  pressure  cases,  a quasi-neutral  fluid  treatment  is  needed. 
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APPENDIX 

GFRED  - A FINITE  GAS  PRESSURE  PARTICLE  CODE 


In  a number  of  IEMP/SGEMP  environments,  the  presence  of  a gas 
can  be  important.  These  environments  include  air  bubbles  in  cables,  residual 

gas  in  systems  at  high  altitudes,  e-beam  cavity  experiments,  to  name  a few. 

19  -3 

In  any  event,  the  number  density  of  molecules  n^  may  range  from  n^  'v  10  cm 
to  essentially  zero.  Finite  gas  pressure  effects  can  cause  the  space  charge 
neutralization  of  an  electron  beam-cloud,  thus  resulting  in  reduced  electric 
fields  and  increased  current  flow.  Or,  at  higher  pressures,  gas  breakdown  is 
possible. 

In  order  to  cover  the  complete  range  of  gas  densities  likely  to  be 
encountered,  it  is  necessary  to  have  a theory  which  describes  phenomena  in  a 
collisionless  plasma  equally  well  as  in  a collision-dominated  plasma  fluid.  It 
is  also  necessary  to  find  a theory  which  goes  smoothly  between  the  two  regimes. 

Before  describing  the  theories  developed,  it  is  worth  considering 
briefly  the  various  phenomena  which  the  theory  should  treat  adequately.  At 
'low'  pressures,  the  beam  electrons  create,  in  passing  through  the  gas,  secon- 
dary electrons.  In  this  regime,  they  accelerate  in  the  space  charge  field  of  the 
beam  electrons  and  escape  to  the  walls  as  runaways.  The  ions  move  much  more 
slowly  and  are  left  to  cancel  some  of  the  space  charge.  For  sufficiently  long 
pulses,  the  space  charge  field  becomes  neutralized.  During  such  a pulse,  the 
electric  field  is  primarily  longitudinal. 

At  higher  pressures,  the  neutralization  takes  place  at  an  early  stage 
in  the  pulse.  Thereafter,  the  electric  field  is  primarily  transverse.  In  this 
regime,  secondary  electrons  will  collide  and  create  further  ionization.  The 
secondaries  will  also  be  scattered  by  elastic  collisions.  As  a result,  it  is 
necessary  to  allow  for  the  possibility  of  secondaries  behaving  either  as  run- 
aways or  as  collision-dominated  drifting  particles.  The  effect  of  the  magnetic 
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field  on  the  electron  motion  (particularly  the  fast  electrons)  may  be  important. 

In  cases  where  the  secondary  electron  density  is  comparable  with  beam  and 
is  characterized  by  a much  lower  temperature,  it  may  be  necessary  to  follow 
effects  such  as  the  two-stream  instability. 

At  yet  higher  gas  densities,  the  ionization  by  secondaries  may  result 
in  avalanche  and  breakdown.  Here,  the  drifting  behavior  will  dominate  and 
the  plasma  behavior  is  fluid-like. 

Proceeding  to  the  highest  gas  densities,  the  collision  rate  for  secondaries 
is  now  so  small  that  they  cannot  gain  enough  energy  between  collisions  to  cause 
ionization.  In  this  regime,  we  expect  the  formation  of  secondaries  to  surround 
the  primary  trajectories.  Recombination  may  now  be  important  and  the  primary 
beam  may,  itself,  be  degraded  and  attenuated. 

To  cover  this  complete  gas  density  range  is  a formidable  task.  It  is 
fair  to  say  that  the  vacuum  case  is  now  well  understood  and  that  the  2-D  codes 
FRED/E,  FRED/EM  and  SQUID  are  all  tested  and  fully  operational.  They 
are  applicable  to  2-D  geometries  of  arbitrary  shap£.  FRED/E  uses  a quasi- 
static model,  whereas  the  other  two  are  electromagnetic  codes.  SQUID  is 
the  most  complete  code,  incorporating  a capability  for  handling  heavily  space 
charge  limited  dynamics  in  an  efficient  manner. 

The  first  steps  toward  a complete  finite  gas  density  capability  have 
been  set  up  with  the  codes  MANFRED  and,  more  recently,  GFRED.  The  former 
is  a finite  density  extension  of  FRED/E,  and  since  it  is  quasi-static,  can  be 
applied  to  very  low  density  cases  where  neutralization  of  the  space  charge  is 
incomplete.  The  more  recent  GFRED  removes  the  quasi-static  limitations  - it 
is  an  extension  of  FRED/EM.  While  applicable  to  a wide  range  of  conditions, 
it  can  not,  in  its  present  form,  be  applied  to  the  entire  gas  density  range  described 
above. 

GFRED  is  a 2-D  axisymmetric,  electromagnetic  particle  code  which 
incorporates,  in  an  approximate  fashion,  a description  of  the  ionization  of  a 
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finite  density  gas.  The  electromagnetic  description  is  carried  out  as  in 
FRED/EM.  The  new  features  are  in  the  treatment  of  secondaries. 

The  method  adopted  to  treat  plasmas  which  may  be  collisionless  or 
collision-dominated  is  to  use  a particle  description  in  which  the  equation  of 
motion  is 


dv 

dt 


e -*■  v -*•  -> 

— (Eh  - x B)  - v\ 
m c 


(A.  1) 


where  v is  the  (energy  dependent)  collision  rate.  The  latter  results  from 

both  elastic  scattering  (at  a rate  v ) and  inelastic  ionization  processes  at  a 

s 

rate  v 


In  GFRED,  as  in  most  particle  simulation  codes,  particles  represent 
large  numbers  of  electrons.  In  the  particle  dynamics,  the  electron  value  for 
e/m  is  used  as  in  Equation  (A.l)  above.  In  calculating  currents  and  charge 
densities,  however,  an  effective  charge  e*  is  used  where  c+/e  is  the  number 
of  electrons  represented  by  that  particle.  The  rate  at  which  a particle  of  charge 
e*  creates  secondary  charge  is  given  by 


^ = n ct  ve*  = v e* 
dt  o I I 


(A.  2) 


where  a is  the  ionization  cross  section.  It  is  assumed  here  that  the  fractional 

ionization  of  the  plasma  is  low  so  that  the  electrons  collide  only  with  neutral 

molecules.  The  collision  cross  section  <7^  is  zero  below  an  energy  threshold 

of  ^20  eV  for  air.  The  energy  dependence  of  a is  based  on  an  analytic  fit 

of  experimental  data  and  is  described  in  Refs.  4,  5 and  6.  The  scattering  rate 

v is  related  to  the  mobility  p by  the  expression 
s 


e 

v - — 
s mp 


and  the  mobility  data  is  provided  in  Ref.  7. 
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For  energetic  particles,  the  equation  of  motion  reduces  to  that  of  a 
particle  in  a collisionless  plasma  because  for  particle  energies  e greater 
than  that  of  the  atomic  ionization  levels  v^\/t  fnt  . On  the  other  hand,  for 
low  energy  particles,  the  inertial  term  is  negligible  and  they  drift  with  a 
velocity 


e 

me 


xB)  = - 


-*■ 

eK 

mr 


(A.  3) 


In  this  limit,  we  have  the  usual  mobility  transport  model. 

This  description  can  thus  be  used  to  go  smoothly  from  collisionless  to 
collision-dominated  regimes  and  to  treat  both  energetic  primaries  and  low-energy 
secondaries.  However,  some  consideration  must  be  given  to  the  high  density 
case  where  the  plasma  is  collision-dominated.  In  this  individual  particle  model, 
secondary  electrons  are  created  on  each  cycle  by  each  of  the  primary  electrons 
in  accordance  with  Equation  (A. 2).  These  particles  are  then  pushed  according 
to  the  requirements  of  Equation  (A.  1).  In  order  to  prevent  the  total  number  of 
code  particles  from  exponentiating  in  this  way,  it  is  necessary,  at  periods 
Tc  = nc  At  (where  n^  is  an  integer  and  At  is  the  time  step)  to  coalesce 
particles  whose  energies  are  such  that  during  the  time  T^  they  do  not  move 
appreciably.  That  is,  they  are  coalesced  if  their  velocity  v is  such  that 


v T <<  Ax 
c 

where  Ax  is  a cell  dimension. 

This  treatment  appears  to  work  well.  Despite  its  simplicity  and  the 
ease  with  which  it  is  incorporated  in  a standard  EM  code  such  as  FRED/EM, 
it  does  contain  most  of  the  physics  necessary  to  describe  the  behavior  of  a 
plasma  before  avalanche  phenomena  are  important.  The  advantages  of  the 
present  model  (over  others)  in  this  area  are  listed  below. 
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1)  'First  principles'  treatment  of  the  secondary  electron 
dynamics  - as  distinct  from  using  a set  of  approximate 
scaling  laws^  to  describe  neutralization,  particle 
drifting,  etc. 

2)  The  theory  of  the  electron  dynamics  goes  smoothly  from 
the  collisionless  to  the  collision-dominated  limits. 

3)  The  model  is  simple  to  incorporate  in  a code  such  as 
FRED/EM  and  calculations  performed  using  it  are  fast. 

Nevertheless,  the  model  has  defects  which  are  evident  at  sufficiently 
high  pressures.  One  of  these  is  the  neglect  of  angular  scattering  during 
collisions.  The  other  involves  the  short  wavelength  electrostatic  oscillations 
which  are  created  when  the  secondary  plasma  density  becomes  very  large. 

Regarding  the  former  of  these,  the  deflection  by  a given  electron- 
neutral  collision  is  not  taken  into  account.  For  calculations  in  which  the 
collision  rate  is  such  that  v At  < 1 for  all  particles,  it  would  be  quite 
straightforward  to  take  account  of  the  details  of  the  collision.  However,  very 
low  energy  particles  may  collide  many  times  during  a time  step.  For  them,  a 
fluid-like  description  is  necessary.  The  coalescing  of  slow  particles  and  the 
treatment  of  their  dynamics  via  Equation  (A.  3)  does,  in  fact,  amount  to  a fluid 
description.  What  is  neglected  is  the  temperature  of  the  fluid.  That  is,  we 
assume  that  the  drift  velocity  v,  which  is  determined  from  Equation  (A. 3)  is 
much  greater  than  the  random  or  thermal  velocity.  While  these  are  good 
assumptions  under  conditions  of  large  fields  (electrostatic  or  inductive)  at  low 
pressures,  they  will  be  invalid  at  sufficiently  high  densities. 

The  other  and  more  important  weakness  of  the  GFRED  treatment 
concerns  the  behavior  of  the  plasma  whenever  the  plasma  period  * 


71 


becomes  smaller  than  the  time-step.  We  saw,  in  Section  3,  how  the  results 
became  unstable  whenever  the  criterion  was  violated.  Use  of  a smaller  time- 
step  is  the  simple,  but  limited,  solution  to  the  problem.  Under  conditions  of 
a high  plasma  frequency,  the  plasma  is  essentially  neutral,  with  space  charge 
irregularities  occurring  on  a length  scale,  v/w  , much  less  than  that  of  the 

r 

system  and  on  a time  scale  much  shorter  than  that  of  other  phenomena  of 
interest.  The  appropriate  model  for  such  instances  is  the  adoption  of  the 
charge  neutrality  condition 

V • J ,0 

to  replace  the  charge  continuity  condition 
If  + V • J = 0 


which  is  inherent  in  GFRED.  Work  is  now  in  progress  to  investigate  such  a 
procedure,  which  when  completed,  will  permit  high  pressure  cases  to  be 
examined. 
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